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ABSTRACT 

Aims. The observation of eclipses during X-rays flares taking place in active cool stars binaries allows us to calculate the position and 
size of the flares. This information cannot be derived by analyzing the decay of the flares, a frequently used approach in the literature 
that requires the assumption of a physical model. We make use of the eclipsing light curve to constrain the set of possible solutions, 
from the geometrical point of view, in two flares of Algol, and one flare in VW Cep. 

Methods. We make use of a technique developed with the system SV Cam (i ~ 90°) and generalize it to binary systems with arbitrary 
inclination. The method simulates all possible geometrical situations that can produce the times of the four contacts of the eclipse. As 
an approximation we assume that the emitting region has a spherical shape that remains unchanged during the eclipse. We however 
show that this is a good approximation for the problem. 

Results. The solutions observed indicate that in two of the three cases the flare cannot be polar (\9\ < 55°) and in a third one the flare 
can be placed either near the pole or at other latitudes. The emitting regions must have a small size (0.002 - 0.5R,), but if interpreted 
as the apex of coronal loops, their length could actually be up to 3. 1 R, for one of the Algol flares. These measurements imply a lower 
limit to the electron density in the emitting region between log n e (cm -3 ) 10.4 and 14.0, and a magnetic field between 70 and 3500 G. 
Similar results are found if the emitting region is assumed to be loop-shaped. 
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1. Introduction 

The determination of the dimension of the X-ray emitting region 
in stellar flares has in most cases been based on assumptions 
on the physics of the flaring emission, in particular the rise and 
decay times. Although it is not possible to directly observe spa- 
tially resolved flares in stars (except for the Sun), it is possible 
to constrain their geometry if the flare is eclipsed by one of the 
stars in the system. 

There are four eclips ed X-ray flares repor ted in the liter- 
ature: one in VW Cep dChoi & Dotanil [19 98). two in Algol 
dSchmitt & Favata|[l999t ISchmitt et al.ll2003l) and very recently 
one in SV Cam dSanz-Forcada et al.ll2006l hereafter Paper I). In 
Paper I we calculated all the possible solutions for the eclipse of 
SV Cam, a system with i ~ 90°, which makes the geometrical 
problem simpler. We extend in this work the methodology ap- 
plied in Paper I to the other three eclipsed flares, which implies 
the use of the equations for the general case i + 90°. 

The first case reported in the literature of an eclipse of a 
flare in X-rays was an event in VW Cep (G5V/K0V), a binary 
system of the W UMa-type. These are eclipsing binaries that 
consist in F-K main sequence stars with short orbital periods. 
They are, in general, strong X-ray sources. The observed X-ray 
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emission is attributed to magnetic activity caused by a dynamo. 
VW Cep is an unusual system, since the primary star (the deep- 
est eclipse, at orbital phase 0=0, occurs when this star is oc- 
culted, with i = 65 °) is the hotter (G5V) but the less massive 
star (M = 0.25 M , R = 0.50 R Q ), while the secon dary (KOV 
M = 0.90 M H , R = 0.93 j? H ) is more massive f see iHilll Il989t 
IChoi & Dotanil 1 19981 and references therein). IChoi & Dotanil 
(1998) analyzed the ASCA light curve of VW Cep proposing that 
the observed dip during the decay of a flare, at <p ~ 0.42 (Fig.[TJ, 
was produced by the eclipse of this flare. From investigation of 
the light curve they also conclude that the flare took place in the 
polar region of the most massive star (the KOV). 

Algol (B8V/K2III) is an interesting target to search for 
eclipses (z = 82.5 °) in X-rays, since the hot component of 
the system can be safely assumed to be dark in X-rays as com- 
pared to the K2III secondary, which is expected to host coro- 
nal features. Two eclipse d flares have been reported in Algol. 
ISchmitt & Favatal (1999) observed, in 1997, an eclipse with 
Beppo-SAX during the decay of a large flare. They estimated that 
the flare reached a maximum of 0.6 R* above the stellar surface, 
in a polar region in order to produce the observed eclipse cen- 
tered arou nd 4> ^ 0.50, and a d ensity of at least ~ 10 11 cm" 3 was 
necessary. ISchmitt et alj d2003l) reported another eclipse during 
a flare in Algol, observed with XMM-Newton in 2001, and they 
used a reconstruction method to conclude that the flare had a 
height of ~ OAR, (0.3 R e ) over the surface (R la = 3.09/3.29/? o ), 
placed near the limb, with a density above ~ 10 11 cm 3 . 

In this work we apply the same technique already used with 
S V Cam (Paper I) to the three cases mentioned above. In the case 
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Fig. 1. VW Cep X-ray light curve (ASCA/SIS and ASCA/GIS) 
in Feb. 2002. The solid lines mark the times of the optical 
eclipses. An eclipse in the flare, around <p ~ 0.45, is displayed 
in the lower panel, where dotted lines indicate the phases of the 
four contacts used in the analysis, with a representation of the 
relative postions of the stars next to them. 



of VW Cep and the XMM observation of Algol we reanalyzed 
the data to calculate the phases of the four contacts. The paper is 
structured as follows: in Sect. 2 we describe the technical details 
of the observations; Sect. 3 develops the geometrical formulae 
used in this technique and the results found for each case. The 
results and their implications in the coronal context are discussed 
in Sect. 4, and the conclusions are given in Sect. 5. 
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Fig. 2. Algol light curve (BeppoSAX/MECS) in Aug 1997. An 
eclipse in the flare, around <p ~ 0.5, is enlarged in the lower 
panel. 



2. Observations 

VW Cep (HD 19 7433) was observed w ith ASCA on 5 Nov. 
1993 23:37 UT dChoi & Dotanil Il998l) . The ASCA satellite 
has four telescopes equipped with two Solid-State Imaging 
Spectrometers (SIS0, SIS1, range 0.4-10 keV), and two Gas 
Imaging Spectrometers (GIS2, GIS3, range 0.8-10 keV). We 
have reanalyzed the light curves by selecting a circular region 
around VW Cep, and subtracting the background level. The light 
curve displayed in Fig. Q]has been obtained by correcting for 
small time shifts between the two GIS datasets and then sum- 
ming the signal. Times corresponding to the four contacts were 
determined in the GIS summed light curve. 

Two observations of Algol (fi Per, HD 19356, HR 936) are 
analyzed in this work. A lgol was observed on 3 Aug. 1997 3:04 
UT with BeppoSAX dSchmitt & Favatalll999l) . We use in this 
case the phases of the four contacts as provided by the authors. 
The corresponding light curve (Fig. |2), was constructed sum- 
ming the data from all the three Medium Energy Concentrator 
Spectrometer (MECS) detectors (range 0. 1-10 keV). 

XMM-Newton observed Algol on 12 Feb. 2002 4:42 UT 
dSchmitt etalJ 120031) . The light curve used for the analysis is 
based on EPIC-PN (range 0.15-15 keV) data that we reanalyzed 
using the standard SAS (Science Analysis Software) version 6.1. 
We determined the four contacts of the eclipse as displayed in 
Fig.0 
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Fig. 3. Algol light curve (EPIC/PN) in Feb 2002. The optical 
eclipse takes place between <p ~ 0.43 and <p ~ 0.57. The eclipse 
in the flare, around </> ~ 0.54, is displayed in the lower panel. 

3. Results 

In Paper I we made an analysis of an eclipse of a flare in 
SV Cam. In that case the high inclination (i ~ 90°) of the sys- 
tem simplifies the equations. Here we calculate new equations 
for the general case i + 90°. We assume that the emitting region 
has a spherical shape, and is defined by the following quantities: 
latitude (9), longitude (A, with origin the meridian of each star 
that is in front of the observer at 0=0), height over the center of 
the star (h) and the radius of the spherical emitting region (Rj,). 
We have made a grid of values for the 4 variables, and we then 
calculated the times when the four contacts of the eclipse take 
place for each set of values. We considered as valid results those 
that agree within 1 time bin (300 s for VW Cep and the Algol 
XMM observation, 600 s for Algol with BeppoSAX) with the 
measured times. We describe here the equations corresponding 
to a system with a primary star that hosts a flare which is eclipsed 
by the secondary. For the simple case (system with i = 90°): 



x'j = —a sin(</>), y\ — 0, z\ — —a cos < 



y' 2 = 0, z' 2 = -z\ 



x' 3 = x\ + h cos(f?) sin(0 + A), y' 3 = y\ +h sin(f5), 
z' 3 =z\ + h cos(f9) cos(0 + A) 

where x', y' correspond to the coordinates in the plane in front 
of the observer (Fig. 0), z' is the coordinate of the axis perpen- 




Fig. 4. Schematic view of the geometry of the stars and the flar- 
ing region (dashed) at phase ~ 0.15. 

dicular to this plane (positive towards the observer), and a is the 
value of the semi-axis of the orbit. x' v yL z\ correspond to the 
position of the emitting region of the flare, and the subindices 1 
and 2 correspond to the positions of the centers of the primary 
and secondary stars respectively. 

Now, for the general case in which i + 90° we need to apply a 
rotation. In the new coordinates system, x,y, z can be calculated 
from the former coordinates: 

x — x' , y = sirn'y' - cos iz', z = cosiy' + sin iz' 

and therefore, 



Z\ — —a cos <p sin i 
Z2 — a cos <p sin i 



xi = — asmq), yi = a cos cos i, 
X2 — a sin <p, yi — —a cos <p cos i, 
xt, — x\ +h cos 8 sin(0 + A), 
y3 = J\ + h[sin6siai - cosf?cos(r/> + A) cos;] 
Z3 =Z\ + /z[sin#cos i + cos#cos(0 + A) sin;'] 

At the first and fourth contacts the distance between the cen- 
ter of the flare and the center of secondary star must be: 



2+R3 )2_ (y3 _ y2) 2 



and similarly, at the second and third contact: 



2 -R i )2-(y 3 -y 2 )2 



Thus, we will have four equations: X3 = Xz + Ai, X3 = 
x 2 + A 2 , X3 = x 2 - A 2 , X3 = x 2 - A t that we can use to cal- 
culate numerically the times of the four contacts for each point 
of the grid. Similar equations can be easily derived for the other 
configurations of the flare. Some geometrical constraints can be 
added to ensure the totality of the eclipse (y$ +R3 < y 2 + R 2 ) and 
a valid position in z (Z3 + Ai < zi)- In cases with large values 
of h we have imposed another constraint: the centrifugal force 
must be lower than the resulting gravitational forces of the two 
stars. In this way, h + R^ must be lower than a certain value that 
is very close to a. Further restrictions can be imposed regarding 
the light curve observed during the phases in which no eclipses 
where observed: (xj, - x 2 ) 2 + (J3 -yz) 2 > Rz +R3 (substitute the 
subindex 2 with 1 to avoid eclipses by the primary). The results 
found are shown in Tables Q] [2] and in Figs.|5]-[Zl More detailed 
results are found in Figs.[8lfTTl and Fig.[T2]show several anima- 
tions of valid solutions in the three flares. 

The physical interpretation of an emitting region shaped as a 
sphere detached from the surface could be the apex of the loops 
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Table 1. VW Cep (G5V/K0V, fl 1;2 =0.50/0.93 Rq) range of possible solutions of the eclipse of a flare in 1993. 



Flaring star 


Eel. star 


en 


AC) 


h (RJ 


h-R, (R«) 


R 3 (Ra) 


logn e (cm -3 ) 


B(G) 


Sec 


Pri 


+ 12.5 - +14.3 


173.3-176.2 


1.53-1.66 


0.65-0.78 


0.002-0.012 


12.8-14.0 


800-3300 


Sec 


Sec 


-17.5 - +49.0 


20.8-25.2 


0.95-1.87 


0.02-1.01 


0.010-0.273 


10.7-12.9 


76-910 






-49.9- -39.1 


20.8-25.2 


1.63-1.87 


0.75-1.01 


0.084-0.223 


10.9-11.5 


88-180 


Sec 


S+P 


-52.7 - -42.9 


20.8-25.2 


1.62-1.87 


0.74-1.01 


0.084-0.197 


11.0-11.5 


98-180 


Pri 


Pri 


+37.5 - +55.0 


20.7-25.2 


0.50-0.54 


0.00-0.08 


0.002-0.014 


12.7-14.0 


710-3000 


All results 


-53 - +55 


20.7-25.2 or 175 


0.50-1.87 
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Fig. 5. Relative positions of the VW Cep stars (dashed lines) during first contact (0=0.337), as seen from the front (left panel) and 
top (right panel). All the possible results are represented with black points indicating the center of the emitting region. A dotted line 
marks the equator (left panel) and the meridian origin of longitude (both panels). Longitude goes positive counter-clockwise. 



involved in the flare, with the emitting region being the place of 
magnetic reconnection. Although the real shape is unlikely to be 
a sphere, such a form is the simplest approximation and would 
be the smallest sphere that contains the real emitting region. We 
have also conducted a test using the shape of a loop (see Sect. 13.41 
and Sect. @}, and another test considering only the central points 
of the eclipse (See Sect.0), yielding similar results. We therefore 
conclude that the use of a sphere gives valid results. 



An additional result can be obtained from the temperature 
(T) and emission measure (EM) of the flare. We can calculate 
the electron density (n e ) of the emitting region by calculating 
the EM during the flare. In this case we use the values of the T 
and EM determined in the original papers when available. Since 
EM~ 0.8 n\ V, it is possible to get the density n e from the net 
flare EM, using V - 4/3nR^. If we consider that the magnetic 
pressure, B 2 /{8n), should be at least as large as the electron pres- 
sure (2n e kT), we can get the minimum magnetic field needed 
to have a stable structure in the flare. In the cases in which the 
whole sphere of the flaring region is not observable, we have cor- 
rected the volume V by a factor accounting for the visible frac- 
tion. The results of the n e and B are listed in Tables [TJ-[2] Next 
we describe the particular details of each of the three eclipses. 



3.1. VWCep 

The eclipse detected in VW Cep with ASCA was analyzed by 
IChoi & Dotanil(ll998l) . They concluded that the dip observed in 
the light curve (see Fig. [TJ corresponds to an eclipse of the flar- 
ing region. Since the decay of the flare lasts almost one orbital 
period, they assume that the flare must occur in a polar region 
because otherwise it w ould be occulted by on e of the stars out of 
the eclipse. Therefore IChoi & Dotanil dl998l) conclude that the 
flare occur in a polar region of the primary. They calculate also 
the size scale of the flare from the ingress and egress times yield- 
ing a size of ~ 5.5 x 10 10 cm (~ 0.8/? Q ), ignoring the rotational 
velocity of the secondary star because the flare is assumed to 
be polar. Rotational modulation will however be produced only 
if the emitting region of the flare is close enough to the stellar 
photosphere to be occulted. 

We have reanalyzed the data, and determined that the phases 
of the four contacts are 0=0.337, 0.3 74, 0.486 and 0.548 respec- 
tively, using the elements provided by Aluig i et aD (1 1 9941) . where 
P orb =0.2783076 and r =2,448,862.5220 (at 0=0 the G5 V start s 
behind the K0V star), and the stellar data given bv lHilll d!989b . 
If we assume that the dip observed corresponds to an eclipse 
(rather than the decay of a flare followed by a second flare), then 
this eclipse partially overlaps with the secondary optical eclipse, 
and therefore it is likely that a flare in the secondary is eclipsed 
by the primary. We have considered also the cases in which the 
flare is self-eclipsed by the host star. If the emitting region is a 
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Table 2. Algol (B8 V/K2III, /?i j2 =3.09/3.29 Rq) range of possible solutions of the eclipses of the flares in 1997 and 2002. 



Observation 


en 


AC) 




h - R 2 (R„) 


Ri (Ro) 


log/j,, (cm J ) 


5(G) 


BeppoSAX 


+ 12.7 - +13.9 


178.6-179.8 


7.81-12.5 


1.37-2.80 


0.033-0.33 


11.8-13.3 


630-3500 




-87.98 - -22.5 


10.10-179.2 


3.14-9.45 


-0.05- + 1.87 


0.26-1.69 


10.8-11.9 


220-750 


BeppoSAX (self-eclipsed) 


+ 1.0- +15.6 


0.4-1.6 


7.50-13.60 


1.28-3.13 


0.29-1.19 


10.9-11.8 


240-690 




-30.0- -16.4 


0.4-1.6 


7.75-13.60 


1.36-3.13 


0.30-1.13 


11.0-11.80 


250-670 


XMM-Newton 


+6.9- +16.0 


197-272 


3.75-7.31 


0.14-1.22 


0.078-0.31 


11.0-11.9 


150-420 




-50.5 - -24.8 


198-323 


5.82-10.7 


0.77-2.25 


0.093-0.76 


10.4-11.8 


76-370 




Fig. 6. Relative positions of the Algol stars (dashed lines) during first contact (0=0.451) in the Beppo/SAX observations, as seen 
from the front (left panel) and top (right panel). All the possible results are represented with black points indicating the center of 
the emitting region. 



flare in the primary separated from the stellar surface, it is not 
possible to distinguish it from a flare in the secondary occulted 
by the primary. The results of the analysis (displayed in Fig. [5] 
and Table [TJ indicate that no flare in polar regions {\9\ > 55°) 
can produce the observed light curve, and that the flare could be 
either associated with the primary or the secondary. Besides, the 
size of the emit ting region is always R -$ < 0.3 R , well below the 
calculations of IChoi & Dotanil dl998l) . 

We have calculated the electron density based on the 
EM duri ng the flare. We hav e used the same values calcu- 
lated bv IChoi & Dolanl (fl998l) : log T u = 6.87,7.35 (K) and 
\0gEM\2 = 52.56,52.58 (cm -3 ). Results are displayed in 
Table [TJ ' 

Although the SIS and GIS light curves are suggestive of an 
eclipse of a flare, we cannot reject the possibility that there are 
actually two flares. However we believe that the slight increase 
in the slope at the "first contact" might be an indication of an 
eclipse rather than a fast decay of a first flare (the SIS light curve 
seemed to indicate more clearly the eclipse interpretation, see 
IChoi & Dotanl[l998h . 

3.2. Algol with BeppoSAX 

BeppoSaX detected a very remarkable flare in Algol in 1997, 
with an increase in X-ray flux by 2 orders of magnitude. Besides, 
the fla re was eclipsed during th e decay in a rather clean light 
curve. Schmitt & F avatal (1 19991) analyzed these data, and they 



concluded that the flaring region was placed near the southern 
pole. If we apply ou r method to the same d ata (using the same 
phases calculated by Schmitt & F avatal fl999l) . we have a set of 
results in a wide range of latitudes, and even some at 8 ~ 15.6° 
in the northern hemisphere. Among the results found, those in 
the northern hemisphere with the flare eclipsed by the primary 
imply a relatively large distance from the secondary star (li), and 
such positions are likely hampered by gravitational forces of the 
secondary star and/or some mass transfer given the low latitude 
and short distance of these solutions to the center of mass of 
the system. However we did not compute the effects of grav- 
ity in order to reject these solutions. Another set of results in 
the southern hemisphere indicates that the flaring region could 
still be present at latitudes between -22 and -59°, in this case 
at lower heights over the secondary, and likely less affected by 
mass transfer processes. Finally, those results related to a self- 
eclipse by the secondary star imply large values of h and low 
latitudes, but no mass transfer can be involved in this case (see 
Table©. 

The phases of the eclipse (0.451, 0.488, 0.505, 0.545) make 
this analysis more complicated than the others. While in the 
other cases the range of solutions does not depend on the ap- 
proximation of a spheric shape for the emitting region, in this 
particular case (totality happens while there is a stellar conjunc- 
tion) it is possible to find other solutions consisting in an eclipse 
of a flare that is very close to the surface. In order to model this 
case we assumed that a sphere partially below the photosphere 
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Fig. 7. Relative positions of the Algol stars (dashed lines) during first contact (0=0.516) in the XMM observations, as seen from the 
front (left panel) and top (right panel). All the possible results are represented with black points indicating the center of the emitting 
region. 



of the star might produce the same light curve, and therefore we 
have relaxed the conditions to allow this geometry to happen. 
This workaround allows us the use of spheric caps in the anal- 
ysis as well. This situation is not considered for the other two 
eclipses (VW Cep and Algol with XMM-Newton) because the 
geometry of the eclipses is not compatible with this configura- 
tion: the Algol case is explained below, and the second contact 
of the coronal eclipse of VW Cep occurs when the photospheric 
eclipse starts (see Fig. and therefore only a region detached 
from the surface can be totally eclipsed at that time. The elec- 
tron density under this configuration has been calculated using 
the volume of the emitting region that is visible (note that in 
some cases the emitting region is partially occulted below the 
secondary star during the whole flare). The solutions found un- 
der this configuration show that the flare could be present at al- 
most any longitude near the south pole of the star, producing the 
same kind of light curve eclipse as the other results mentioned 
in former paragraph. 

3.3. Algol with XMM-Newton 

XMM-Ne wton detected a flare in Algol in 2002 that was ana- 
lyzed by Sch mitt et al.l (120031) . who concluded that an eclipse 
of the flare was taking place, indicating that the flare should 
have occurred near the limb in the northern hemis phere since 
the fir st contact takes place at (p >0.5. In this case Schmitt et al. 
(2003) used a reconstruction of the flaring plasma region using 
the whole light curve, concluding that the flare takes place at a 
height of ~ 0.1/?* (R 2 = 3.29fl G ) and 9 < 69°. 

We have determined the phases of the four contacts at 0.516, 
0.527, 0.554, 0.565 using the ephemeris of |A1-Naimiv et al.l 
(fl985l) . where P orb =2.8673285 d and r =2,445,7 39.003. These 
values are in agreement with those proposed by Schmit fet alj 
(2003), although they only expressed the first three contacts 
and assumed that the fourth contact takes place at the end of 
the observation. We do not expect this to affect the range of 
solutions. However there is an alternative interpretation of the 
light curve: the eclipse could actually be starting at phase 0.506 



(the apparent peak of the flare), while the change in slope at 
4> ~ 0.5 16 could actually correspond to the flare maximum emis- 
sion. The determination of the fourth contact has also some un- 
certain ty. The analysis of the light curves of different energy 
bands (Schmitt et al. 2003) seems to suggest that the eclipse fol- 
lowing a flare maximum is the best option, so we will refer to 
this situation from now on. 

Given the phases involved in the eclipse (see Fig.0 it is eas- 
ier to analyze this flare than the BeppoSAX flare, since a config- 
uration of an emitting region overlapping with the star is not pos- 
sible: the eclipse starts after <p ~ 0.5, invalidating the solutions in 
the southern hemisphere, and those in the northern hemisphere 
were rejected in the analysis. We only found valid solutions for 
the configuration in which the flare (hosted by the secondary) is 
eclipsed by the primary star. Solutions, including the result re- 
ported by Sch mitt et al.l d2003l) . are possible both in the northern 
(7° < 9 % 16°) and southern hemisphere (-50° < 9 < -25°), 
and given the position of the flare (Figs. [7] [TTJ in most cases it 
is unlikely that the flaring region is affected by the mass transfer 
between the two stars. 

Since lSchmitt et al.l(l2003l) do not report the value of the EM, 
we made a spectral fit to the spectrum in the quiescent emission 
(log EM [crrT 3 ] =53.64) and in the flaring non-eclipsed emission 
(log EM [ctrT 3 ]=53.89), resulting in a net value for the flare of 
log EM [ctrT 3 ]=53.53 at log T [K]= 7.50. The range of possi- 
ble values of the electron density (Tab le O are c onsistent with 
those reported bv lSchmitt et al.l d2003l) . but unlike ISchmitt et al.l 
(2003) we also found solutions in the southern hemisphere, and 
a wider range of solutions in the northern hemisphere. 

3.4. Eclipses of coronal loops 

The shape chosen for the analysis has been the simplest possi- 
ble, a sphere. The use of a different geometry will further re- 
strict the number of solutions. Since there is an infinite range 
of shapes possible, we have chosen only one particular case to 
show the deviation from the solutions found with the sphere. We 
chose a shape similar to what is observed in solar coronal loops. 
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Fig. 13. Schematic view of the coronal loop. The loop lies around the XZ plane. The Z axis points opposite from the center of the 
host star and Y axis is paralel to the closest meridian. 



Following the notation used in the paper we have defined a torus 
contained in a sphere of radius R3 (the external radius of the 
torus), a cross-section radius fixed to R4 = 0.2 /?3, and a circular 
curvature. The arc is placed perpendicular to the stellar surface 
(see Fig.fOl). with a local reference frame with the Z axis point- 
ing opposite from the center of the star and the Y axis parallel 
to the tangent of the nearest meridian. In this local reference 
frame we can define each of the points (xuyuZi) of the loop us- 
ing two angles (B, 5) and a radial variable (r, the radius of the 
cross-section) as follows: 

x\ = (Rj, - R4 + r cos 6) cos B, y\ = r sin 6, 

Zi = (Rj - Ra + r cos 5) sin B 

where r varies from to R4, B from 0° to 180° and 6 from 0° 
to 360°. In order to account for the eclipses coverage we need to 
transform the local coordinates (xuyuZi) into the coordinates of 
the rest reference frame used for the stellar system (Fig. [4]). We 
need to apply 3 rotations, in latitude (9), in longitude (A + <p), and 
in inclination (/): 

a rotation in latitude (around X), 

x = x, y = (cos 6)y + (sin 8)z, z' = (— sin ff)y + (cos 8)z, 

a rotation in longitude (around Y), 
x" - cos(0 + A)x' + sin(0 + A)z', y" - y', 
z" = - sin(</> + A)x + cos((p + A)z, 

and a rotation due to inclination (around X), 
x" = x" , y'" = (sin i)y" - (cos i)z" , 
z" - (cos i)y" + (sin i)z", 



Therefore, the position of a point of the loop in the same 
reference frame used in the problem will be: 

x[" = JC3 + cos(0 + A)x\ + sin(</> + A) (- sin 9yi + cos 0z\), 

y'i" - + sm ' ( cos 6)>i + sin Ozi) - cos i (- sin(0 + A) %\ + 

cos(0 + A) (- sin 6y\ + cos 0z\), 
z[" = Zi + cos i (cos Oyi + sin 6zi) + sin i (- sm(<p + A) Xi + 

cos(0 + A) (- sin 6y\ + cos 0z\), 



We have applied this method to search for the results in the 
case of the flare observed with XMM-Newton in Algol. In Fig. [14] 
we display a valid solution near the equator, with a simulation in- 
cluded in Fig. [12] The possible solutions (Fig. [15]) are very sim- 
ilar to those found for a spheric shape (the h and Rj, are slightly 
different since R3 has a different definition). The density derived 
from this geometry would also be higher than that derived from 
the spheric case. 



4. Discussion 

The approach presented here constrains the geometry (location 
and size) of the flaring region for eclipsed flares. A simple calcu- 
lation also yields a lower limit of the electron density (the plasma 
does not need to be homogeneously distributed in a sphere). 
The results are consistent with the range of values derived from 
density-sensitive line ratios in active stars. The method assumes 
that the emission measure remains constant during the eclipse 
(this is important only for the calculation of electron density and 
magnetic field), and that the spherical shape of the emitting re- 
gion remains unchanged during the eclipse. This shape has been 
chosen as the simplest possible approximating a compact region. 
Spheres with a filling factor lower than 1 would imply higher 
values of the electron density, and we do not expect that other 
shapes vary substantially the results of latitude, longitude and 
height, except in cases like the BeppoSAX Algol flare, where the 
flaring region can be actually attached to the surface and have a 
significant surface filling factor. The solutions where the emit- 
ting region is detached from the surface could correspond to the 
apex of the loops involved in the flare, although magnetic re- 
connection in the solar flares does not happen necessarily at the 
apex, and the site of reconnection does not need to be the site of 
maximum X-ray emissivity. Although it is unlikely that the emit- 
ting region is a sphere, we consider this a good approximation 
in order to search for the location. We have conducted two tests 
to quantify the relevance of the shape chosen with respect to the 
real solutions. A first test was conducted using a coronal loop 
instead of a sphere (see Sect. 13.41 ), and considering the four con- 
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Fig. 14. An example of a valid solution using a coronal loop placed at 9 — 10.1°, A = 258.2°, = 0.53 R and h = 3.52 R . Views 
from the front (left panel) and the top (right panel) during first contact. 



tacts of the eclipse. The orientation of the loop (with the shape of 
a torus) was fixed, as well as the cross-section radius (0.2 ^3). 
The set of results is similar to that of spherical shape, but the 
range of solutions is smaller. In the second test, we have con- 
sidered the simpler case in which only two contacts are used in 
the model. These contacts would correspond to the ingress (mid- 
point between <p\ and #2) and egress (midpoint between (pj and 
04 ) of the center of the emitting region, and therefore indepen- 
dent of its shape. The test conducted has shown that the range of 
solutions is similar to those obtained for the restrictions imposed 
by the use of four contacts and a spherical shape. The results are 
displayed in Figs. fT6HT8l These two tests confirms that the ap- 
proximation of a sphere can be valid also as an approximation 
for more complicated shapes. 



While th e original analyses of these events pointed to wards 
polar flares, dChoi & Dotanilll998tlSchmitt & Fav ata 1999), this 
work has shown that other configurations are possible from the 
geometric point of view. A polar location is only possible for 
the Algol flare observed with BeppoSAX. In this work we do 
not privilege any solution among the possible results, and con- 
sider all configurations as equally possible as long as they do 
not violate any physical law. Gravity forces or mass transfer be- 
tween close binaries should be considered in order to constrain 
the number of solutions. 



5. Conclusions 

We have systematically searched the solution space to calcu- 
late the geometrical characteristics of the emitting region of four 
flares in active stars. The method uses the times of the four con- 
tacts of the eclipse and simulates all possible geometrical situ- 
ations assuming that the emitting region has a spherical shape 
and no substantial geometrical changes during the eclipse. Two 
tests (including the use of a loop) were conducted to prove that 
the approximation of a spherical shape is appropriate, even if the 
real shape is not a sphere. The solutions found in the three flares 
analyzed in this work and that of Paper I show that polar flares 
are not possible in three of the four flares, and electron densities 
must be larger than log n e (cm _3 )=10.4, consistent with measure- 
ments from line ratios of density sensitive lines in active stars. 
The magnetic fields needed to confine these loops range in both 
systems between 75 G and 3300 G. The emitting regions have 
sizes that range from 0.002 to 0.5 R*, and they can be at a dis- 
tance (interpreted as loop length) of up to 3.1 R, from the star's 
surface. Further refinement based on other physical limitations, 
such as gravity or mass transfer could further constrain the set of 
solutions. 

Acknowledgements. This research is based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and contributions directly 
funded by ESA Member States and NASA. This research has made use of 
NASA's Astrophysics Data System Abstract Service. JS acknowledges the sup- 
port by the ESA Research Fellowship Program. 



The analysis poses a clear lower limit in electron density of 
log n e (cm _3 )= 10.4, well in agreement with the minimum val- 
ues determined from line ratios analysis in active stars. The 
largest values compatible with the geometrical configurations, 
specially those with log n e (cm -3 ) k, 13, are larger than those 
commonly observed in th e X-ray line ratios of most stars (e.g . 
ISanz-Forcada et alj 120031; iTesta et alj 12004 iNess etalJ |2004 . 
The measurements of the size of the emitting region shows that 
although the emitting regions are small compared to the stellar 
size, we cannot reject large loop lengths (measured through the 
variable h) by applying only geometrical constrains based on the 
eclipse contacts. 
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Fig. 8. Results of the flare in VW Cep, flare in the secondary. Left: eclipsed by primary. Right: self-eclipsed by secondary. 
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Fig. 9. Results of the flare in VW Cep. Left: flare in primary, self-eclipsed by primary. Right: flare in secondary, the eclipse starts 
behind secondary and ends behind the primary. 
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Fig. 10. Results of the flare in Algol observed with BeppoSAX, solutions in the southern hemisphere. Left: The emitting region is 
detached from the surface of the star (the results around A ~ 1° correspond to self-eclipse by the secondary). Right: The emitting 
region can overlap with the star (therefore it can be a spheric cap). 
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Fig. 11. Results of the flare in Algol. Left: Observed with BeppoSAX, solutions in the northern hemisphere (the results around 
A ~ 1° correspond to self-eclipse by the secondary). Right: Observed with XMM-Newton. 



Fig. 12. Several examples of animations of valid results with different configurations. A white circle (or loop) indicates the emitting 
region of the flare being eclipsed. The latitude and longitude used in the simulations are shown in the upper panel, while the light 
curve is displayed in the lower panel. 
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Fig. 17. Results for the case of two contacts only in VW Cep. From top to bottom: flare in secondary eclipsed by primary, self- 
eclipsed by secondary only, and eclipse starting behind secondary and ending behind primary. 
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Fig. 18. Results for the case of two contacts only. Upper panels: VW Cep, self-eclipsed by primary. Lower panels: eclipse in Algol 
observed by XMM (flare in secondary, eclipsed by primary). 



